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Generalized statistical mechanics for superstatistical systems 

Christian Beck 
School of Mathematical Sciences, Queen Mary University of London, Mile End Road, London El 4NS, UK 

Mesoscopic systems in a slowly fluctuating environment are often well described by superstatis- 
tical models. We develop a generalized statistical mechanics formalism for superstatistical systems, 
by mapping the superstatistical complex system onto a system of ordinary statistical mechanics with 
modified energy levels. We also briefly review recent examples of applications of the superstatis- 
tics concept for three very different subject areas, namely train delay statistics, turbulent tracer 
dynamics, and cancer survival statistics. 



I. INTRODUCTION 



o 

Nonlinear dynamical processes often create a fluctuating environment for a given mesoscopic system jl'] . This leads 
\Q [ to mixing of the mesoscopic dynamics and that of the environment. If there is sufficient time scale separation, then 
very often superstatistical models yield a good effective description. The superstatistics concept has established itself 
as a powerful tool to describe quite general classes of complex systems [H-llll [23 - r26J | . The basic idea is to characterize 
the complex system under consideration by a superposition of several statistics on different time scales, for example 
one corresponding to ordinary statistical mechanics (on a mesoscopic level modelled by a Langevin equation) and the 
other one corresponding to a slowly varying inverse temperature field f3(x,t) or some other relevant parameter. 

There may be either spatial or temporal variations of the environment. The environment is represented by a suitable 
parameter entering the stochastic differential equation describing the mesoscopic system. The superstatistics concept 
can be applied in quite a general wa y, a nd a couple of interesting applications for a variety of complex systems have 
been pointed out recently [12h21L I27.. 28]. Essential for this approach is the existence of sufficient time scale separation 
so that the system has enough time to relax to a local equilibrium state and stay within it for some time. 

The stationary distributions of superstatistical systems, obtained by averaging over all /3, typically exhibit non- 
Gaussian behavior with fat tails, which can be a power law, or a stretched exponential, or other functional forms as 
well Q. In general, the superstatistical parameter /3 need not to be an inverse temperature but can be an effective 
parameter in a stochastic differential equation, a volatility in finance, or just a local variance parameter extracted fronr 
some experimental time series. There are interesting applications in hydrodynamic turbulence [21 20l I29ll30| . for defect 
turbulence [12l |. for cosmic rays 11311 a nd other scattering processes in high energy physics I3ll l32j [. solar flares (14J . 
share price fluctuations [la . 1271 l33l l34| , random matrix theory |16l l28l . l35| , random networ ks l36l| , multiplicative- noise 



stochastic processes [37( , wind velocity fluctuations [17|, [18( , hydro-climatic fluctuations [19| , the statistics of train 
departure delays |38l ] and survival statistics of cancer patients [39] . Maximum entropy principles can be generalized 
in a suitable way to yield the relevant probability distributions that characterize the various important universality 

q^ [ classes in superstatistics 0, Hfl-Ej] . 

In this paper we shall develop a new theoretical approach to superstatistics, by formally mapping the superstatis- 

i^* tical system onto a system of ordinary statistical mechanics where the energy levels are modified in a suitable way. 

£— n i This approach yields a new interesting theoretical tool to further develop the generalized statistical mechanics of 
superstatistical complex systems, and is described in detail in section 3. We also briefly review some recent examples 

t— I . of applications of superstatistical techniques. Our three examples, all from very different subject areas, are train 
delays on the British railway network, velocity signals in hydrodynamic turbulence, and the survival statistics of 
cancer patients. 



II. REMINDER: WHAT IS SUPERSTATISTICS? 



The concept is best illustrated by starting with a particular example of superstatistics, in fact the one that was 
considered first in [lCt lll| . Consider the following well-known formula: 
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FIG. 1: A spatially inhomogeneous situation of mesoscopic systems (sketched as circles) embedded into a fluctuation environ- 
ment with different inverse temperatures fy. A Brownian particle moves through the different regions with different inverse 
temperatures. 

is the x 2 ( or T) probability distribution and fto and q are parameters (q > 1). 

We see that averaged ordinary Boltzmann factors e~@ E with x 2 -distributed /3 yield an effective Boltzmann factor of 
g-exponential form, given by the right-hand side of eq. ([T]). The physical interpretation is that nonequilibrium systems 
with temperature fluctuations give rise to an effective description in terms of more general Boltzmann factors. In 

statistics 



the x -distribution was advocated for f(/3), because at that time the aim was to better understand q- 
44 1 from a dynamical point of view. General /(/3) were then suggested in [l|. In that paper also the name 



'superstatistics' was created. This name was simply an abbreviation for the fact that there is a superposition of two 
(or several) statistics. In no way this name wants to indicate that this type of statistics is 'superior' to others. 

One can also construct dynamical realizations of superstatistics in terms of Langevin equations with parameters that 
fluctuate on large time scales [Tljj . These local Langevin equations decribe the mesoscopic system under consideration. 
The situation is sketched in Fig. 1. The simplest example would be locally a linear Langevin equation 



v = — 7i> + aL(t) 



(3) 



with slowly fluctuating parameters 7, a. Here L(i) denotes Gaussian white noise. This describes the velocity v of a 
Brownian particle that moves through spatial 'cells' with different local (3 :— j/(2a 2 ) in each cell (a nonequilibrium 
situation). If some probability distribution f(/3) of the inverse temperature /3 for the various cells is given, then 
the conditional probability given some fixed (3 in a given cell is Gaussian, p(v\f3) ~ e~^ v , the joint probability is 
p(v,P) = f((3)p(v\(3) and the marginal probability is p(v) = / °° f((3)p(v\(3)d(3. Integration over (3 yields effectively 
Boltzmann factors that are more general than Gaussian distributions, which depend on the specific properties of /(/3). 
If there are only finitely many cells, then the integral is understood to approximate the average over a large number 
of cells. 

The principal idea of superstatistics is to generalize this example to much broader systems. For example, j3 need 
not be an inverse temperature but can in principle be any intensive parameter. Most importantly, one can generalize 
to general probability densities f(j3) and general Hamiltonians. In all cases one obtains a superposition of two different 
statistics: that of /3 and that of ordinary statistical mechanics. Superstatistics hence describes complex nonequilibrium 
systems with spatio-temporal fluctuations of an intensive parameter on a large scale. The effective Boltzmann factors 
B(E) for such systems are given by 



B(E) 



f(/3)e-^df3 



(4) 



A lot of research has been done in this direction in recent years. If there is locally Gaussian behaviour, then 
the theory of superstatistics is clearly related to the theory of Gaussian scale mixtures. More generally one can 
prove a superstatistical generalization of fluctuation theorems [3( , develop a variational principle for the large-energy 



asymptotics of general superstatistics B, proceed to generalized entropies for general superstatistics [a, EJ |43[ , let 
the q- values in eq. ((TJ fluctuate as well [7[ , and prove superstatistical versions of a Central Limit Theorem Q . There 
are also relations with fractional reaction equations [43], random matrix theory [la, [28|,|35J], networks |3a]. and path 
integrals [f| . Very useful for practical applications is a superstatistical approach to time series analysis [2, [is ^7] ■ 
Applications have been poin ted out for 3d hydrodynamic turbulence [2], [20J, [29J, [30j , wind velocity fluctuations [17J, [18| , 
finance and economics [2l|, |3J, 0, 1481. b linking quantum dots (49j, cosmic ray statistics [l3| and quite generally 
scattering processes in particle physics [3ll I32J . The concept has also been useful to analyze hydroclimatic fluctuations 
19f as well as the statistics of train delays on the British railway network [38J . There are also medical applications 



III. MAPPING SUPERSTATISTICS ONTO CONVENTIONAL STATISTICAL MECHANICS 

Consider a system of ordinary statistical mechanics with energy levels Ei of microstate i. We are looking at a 
canonical ensemble and a priori the inverse temperature /3 is fixed. Now look at identical copies of the system but 
with different temperatures /3j in each spatial cell j, at a given snapshot of time. This is a nonequilibrium situation. 

Let /3o = L f(f3)(3d/3 be the average inverse temperature. We may formally consider a super-Hamiltonian describ- 
ing the entire system which in the different spatial cells has effective energy levels Ei, by writing 

l3 Q E t =l3,E t . (5) 

Apparently this means the super-Hamiltonian has energy levels Ef given by 

Ef = if-* (6) 

Po 

in cell j. 

Since ordinary statistical mechanics is valid for arbitrary energy levels, in particular also for the E\ , we may now 
do ordinary statistical mechanics for the super-Hamiltonian and introduce the partition function Z((3q) of the entire 
system as 



Z(M = 5>-A*i" (7) 
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f(P)e-P"'dl3, (9) 

where in the last step the sum over j is approximated by an integral. Since Z(/3o) — X)i i e ~ i is an ordinary 
partition function (though with exotic, locally modified, energy levels), it is now possible to do ordinary statistical 
mechanics for this superstatistical nonequilibrium system, with all the known formulas. 

We regard the free energy of the superstatistical system as a function of the mean inverse temperature po and define 
it as 

F(/3 ) = -±-logZ(/3 ). (10) 

Po 

In the statistical mechanics formalism it is often convenient to work with the function ^(^o) := (3 F(^ ). Defining 
kTn = 4- one has 

F = U-T S (11) 

* = [3 U-S (12) 

But one has to be careful here what the meaning of the symbols U and S is: U is now the mean energy of the 
energy levels Ei, rather than Ei, and indeed this means that U is a global mean energy corresponding to the entire 
superstatistical system consisting of many cells. One has 



U = 5> W 4 W (13) 



« 4- r wmpuiooiP), (15) 

Po Jo 

where Ui oc (/3) = ^2piEi is the local internal energy in a cell of inverse temperature /?. The entropy is still given by 
the Boltzmann-Gibbs-Shannon form, but formed with the exotic energy levels E\ : 



S = -k^p^logp^ (16) 



= - fe EwT rA? (- lo ^(«-Af) ( 17 ) 

it (A}) 

= f3 U + log Z{0 O ) (18) 

In this way we have formally mapped the superstatistical nonequilibrium system onto an (exotic) equilibrium system 
of ordinary statistical mechanics with average inverse temperature po and a new type of Hamiltonian, corresponding to 
the energy levels E\ . We should remark that the above idea of mapping superstatistics onto the statistical mechanics 
of an exotic Hamiltonian is completely new and different from previous attempts of developing a generalized statistical 
mechanics for superstatistical systems [3 [23, So, EH • 

IV. POSSIBLE SUPERSTATISTICAL DISTRIBUTIONS /(/?) 

The distribution /(/3) is determined by the dynamical large-scale structure of the complex system under consider- 
ation. There have been attempts to derive the specific form of /(/3) relevant for a given complex system wi th g iven 
constraints from a generalized maximum entropy principle. We don't elucidate this further here but refer to |43l | and 
references therein for further details. Actually, what we want to do here is to proceed to practical applications. The 
relevant question is what type of /(/3) are typically seen for experimental data as generated by a generic complex 
system. 

There seem to be three different superstatistics that are of utmost importance Q . These are (a) x 2 -superstatistics 
(= Tsallis statistics), (b) inverse x 2 -superstatistics, and (c) lognormal superstatistics. 

In case (a), /(/3) is given by the Gamma distribution 

1 / \ ™/ 2 



m = WT){wJ "' e ~ m - (1<J) 

where again /3q is the average of f3. This generates generalized Boltzmann factors B(E) that decay with a power law. 
n is a parameter characterizing the number of degrees of freedom. 
In case (b), /(/3) is given by 



/(« = ^m p- /2 - 2 e-^. (20) 



n / n \ Till 

r(§) 

In this case the generalized Boltzmann factors B(E) ~ j f((3)e~P E decay as e~@^ E for large E. 
Finally, in case (c) /(/3) is given by the lognormal distribution 

1 f-(ln^) 2 l 

/( «^^5 exp {^^r (21) 

where /i and s are suitable parameters. In the remaining sections, we briefly describe one example for each of these 
three different cases. 

V. TRAIN DEPARTURE DELAYS 

Traffic delays on the British railway network are reasonably well described by x 2 -superstatistics. The probability 
density of observed train departure delays of length t has been analyzed in detail in [381 ] . Millions of departure times 
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FIG. 2: Observed departure delay statistics on the British railway network (data from [3g|). The solid line is a q-exponential 
fit as given by eq. ()23[) . 



were automatically stored and evaluated. The Oth-order theoretical model for the waiting time t is a Poisson process 
which predicts that the waiting time distribution until the train finally departs is P(t\/3) — /3e - ^*, where f3 is some 
parameter. But this does not agree with the actually observed data 38]. A much better fit is given by a g-exponcntial, 
see Fig. 2. 

What may cause this power law that fits the data? The idea is that there are fluctuations in the parameter j3 as 
well. These fluctuations describe large-scale temporal or spatial variations of the British rail network environment, 
which take place on a much larger time scale than the actual train departures. /3-fluctuations are e.g. produced at the 
begin of the holiday season with lots of passengers, or if there are problems with the track or bad weather conditions. 
Also there can be extreme events such as derailments, industrial action, terror alerts, etc. The observed long-term 
distribution of train delays is then a mixture of exponential distributions where the parameter /3 fluctuates: 



f(f3)p(t\p)dp = / /(/3)/3e-^. 



Pit) 



For a x 2 -distributed /3 with n degrees of freedom one obtains 

p(t) = C-(l + b(q-l)t)^ 

where q = 1 + 2/(n + 2), b = 2/3o/(2 — q) and C is a normalization constant. 



(22) 



(23) 



VI. TURBULENT FLOWS 



Various aspects of hydrodynamic turbulence are quite well described by lognormal superstatistics. In this case 
the mesoscopic local dynamics corresponds to a single tracer particle that is advected by the turbulent flow. The 
environment of the tracer particle changes. For a while it will see regions of strong turbulent activity, then it will move 
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FIG. 3: Measured probability density of small-scale velocity differences in Lagrangian turbulence (data from [50]) and a fit of 
the form (I25p (see [3(| for more details). 



on to calmer regions, just to continue in yet another region of strong activity, and so on. This is a superstatistical 
dynamics similar to Fig. 1. In 'Lagrangian turbulence' one is interested in the statistics of velocity differences 
u(t) := v(t + t) — v(t) of the particle on a small time scale r. For r — > this velocity difference becomes the local 
acceleration a(t) = u(t)/r. A superstatistical Lagrangian model for 3-d velocity differences of the tracer particle has 
been developed in [3CJ ■ The mesoscopic dynamics is a superstatistical Langevin equation of the form 



-ju + Bri x u + oL(t) . 



(24) 



Here 7 and B are constants. The term proportional to B introduces some rotational movement of the particle, 
simulating vortices in the flow. The noise strength a and the unit vector n evolve stochastically on a large time scale 
T a and T r -j, respectively. T a is of the same order of magnitude as the integral time scale Tl, whereas 7 -1 is of the 
same order of magnitude as the Kolmogorov time scale r n . In this model the Reynolds number R\ is basically given 
by the time scale ratio T a ^ ~ Tl/t^ ~ R\ » 1. The time scale Tft >> r n describes the average life time of a region 
of given vorticity surrounding the test particle. Further details are described in [301 ] . 

The parameter /3 is again defined to be /3 := 27/er 2 , but it does not have the meaning of a physical inverse 
temperature in the flow. Rather, one has /3 _1 ~ i/ 1 ' 2 (e) _1 ' 2 e, where v is the kinematic viscosity and (e) is the average 
energy dissipation, e is known to fluctuate in turbulent flows. Kolmogorov's theory of 1961 suggests a lognormal 
distribution for e, which leads us to lognormal superstatistics. For very small r the 1-d acceleration component of the 
particle is given by a x = u x /t and one gets out of the model the 1-point distribution 
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(25) 



which agrees very well with experimentally measured data of the acceleration statistics (Fig. 3). 

The 3-d superstatistical model of [30] predicts correlations between the three acceleration components. An intrinsic 
property of the model is that the acceleration a x in x direction is not statistically independent of the acceleration a y 
in y-dircction. We may study the ratio R := p(a x ,a y )/(p(a x )p(a y )) of the joint probability p(a x ,a y ) to the product 
of 1-point probabilities p(a x ) and p{a y ). For independent acceleration components this ratio would always be given 
by R = 1, whereas the 3-d superstatistical model yields the prediction 



R 



j °° pme-^^+^dp 



r °° /3i/2/(/3)e-^ 2 ^/3 / °° /3V2/(/3)e 



5d/3 



(26) 



The trivial result R = 1 is obtained only for /(/?) = 6(0— Pq), i.e. no fluctuations in the parameter j3 at all. Fig. 4 shows 
R := p(a x ,ay)/(p(a x )p(a y )) as predicted by lognormal superstatistics: Experimental measurements of acceleration 
correlations yield very similar results to those predicted above [30, [51j . 
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FIG. 4: Predicted shape of correlations between acceleration components in a turbulent flow as given by eq. (|26[) . /(/?) is the 
lognormal distribution. 



VII. SURVIVAL STATISTICS OF CANCER PATIENTS 

Data of the survival statistics of cancer patients can be well fitted using models based on inverse x 2 -superstatistics. 

A superstatistical model of the progression of metastasis and the corresponding survival statistics of cancer patients 
has been developed in [39J ■ The final result that comes out of that model is the following prediction for the probability 
density function of survival time t of a randomly chosen patient that is diagnosed with cancer at t = 0: 



P(t) = 
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This can also be written as 
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(28) 



where K v (z) is the modified Bcssel function. Note that this is inverse x 2 -superstatistics. The role of the parameter /3 
is now played by the parameter A, which in a sense describes how aggressively the cancer propagates. This parameter 
has fluctuations from patient to patient. 

The above formula based on inverse x 2 -superstatistics is in good agreement with real data of the survival statistics 
of breast cancer patients in the US. The superstatistical formula fits the observed distribution very well, both in a 
linear and logarithmic plot (see Fig. 5). 

When looking at the time scales in the above figures one should keep in mind that the data shown are survival 
distributions conditioned on the fact that death occurs from, cancer. Many patients, in particular if they are diagnosed 
at an early stage and treated accordingly, will live a long healthy life and will die from something else than cancer. 
These cases are not included in the data shown above. 



VIII. CONCLUSION AND OUTLOOK 



In this paper we have dealt with mesoscopic and other complex systems that are embedded into a temporally 
changing or spatially fluctuating environment. If there is sufficient time scale separation, then a mixture of differ- 
ent statistics (a superstatistical description) is an appropriate method to describe these types of complex systems. 
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FIG. 5: Survival time statistics of breast cancer patients diagnosed with cancer at t — 0, both in a linear and double logarithmic 
plot. The solid line is the superstatistical model prediction [39l |. 



Basically this means that one mixes ordinary statistical mechanics with another statistics of e.g. local temperature 
fluctuations. 

In the first part of the paper we have pointed out how superstatistical complex systems can be mapped onto systems 
of ordinary statistical mechanics. The key point is that one deformes the effective energy levels in a suitable way and 
then applies the well-known techniques of ordinary statistical mechanics to this (exotic type of) super-Hamiltonian. 

In the second part of the paper we summarized a few recent applications of the superstatistical approach to real- 
world problems, which covered quite a range of different subject areas. We studied train delay statistics, turbulent 
tracer dynamics, and survival statistics of cancer patients. Many other areas may benefit from a generalized statistical 
mechanics formalism for superstatistical systems as well. 

In a year where there is a good chance to finally experimentally confirm the long-awaited Higgs particle, it might be 
appropriate to end this paper by mentioning that scattering processes in high energy physics can also be well-described 
by superstatistical models. The experimentally observed power laws of differential cross sections and energy spectra 
at very high energies have been modelled in terms of superstatistical generalized statistical mechanics [Tj, l3lll32l . I52J . 
Superstatistical techniques have also been recently used to describe the space-time foam in string theory [53j, and a 
generalized statistical mechanics model underlying chaotic types of vacuum fluctuations yields a Higgs mass prediction 
of 154 GeV [5J, |55j. It seems there is a lot of scope for relevant contributions of generalized statistical mechanics in 
high energy physics and quantum field theory. 
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